########Exercicio 6.1#############
###Crie seus dados
conj1<-rnorm(10, mean=6, sd=3)
conj2<-rnorm(10, mean=7.5, sd=3.2)
source("simula.r")
dif<-mean(conj1)-mean(conj2)

#a diferenas absoluta entre as mdias  1.311, h uma diferena 
#entre as mdia, mas no sabemos se essa diferena  gerada pelo acaso.

#gerando uma simulaao

sim.dif<-simula(conj1,conj2, nsim=1000)
dados.maior<-sum(sim.dif$diferencas>=dif)
dados.menor<-sum(sim.dif$diferencas<=dif)
pvalue<-(dados.maior/length(sim.dif$diferencas))+
  (dados.menor/length(sim.dif$diferencas))
#as amostras no so diferentes

#teste de mdias
sim.maior<-simula(conj2,conj1,nsim=1000,"uni")
pvalue2<-(sum(sim.maior$diferencas>=dif))/length(sim.maior$diferencas)
#testt
t.dif<-t.test(conj1,conj2,alternative="two.sided")
t.maior<-t.test(conj2,conj1,alternative="greater")

#no existe diferenas entre as amostras

####Exerccio 6.2#########

mudas=read.table("altura-mudas.csv", header=T, sep=",")
str(mudas)
head(mudas)
tail(mudas)
mudas$substrato
mudas$especie
mudas$bloco
mudas$altura
#1.
#a varivel medida  a altura das especies. Como fator temos as especies
# e o substrato.
#temos duas rplicas do fator especie e dez rplicas do fator substratos
#transformando substrato em fator
substrato<-as.factor(mudas$substrato)
mudas$substrato<-substrato
str(mudas)

#ANOVA na unha
tamboril<-mudas[mudas$especie=="tamboril",]
tamboril
#somatrio dos desvios quadraticos geral dos grficos
media.altura=mean(tamboril$altura)
dif.geral<-tamboril$altura-media.altura
ss.total<-sum(dif.geral^2)
#media de altura de cada substrato
media.substrato<-c(tapply(tamboril$altura,INDEX=tamboril$substrato, 
                                FUN=mean))
#Somatrio do desvio quadrtico intragrupo
str(tamboril)
head(tamboril)
tamboril
vetor.media.substrato<-rep(media.substrato,each=6)
tamboril$media.grupo<-vetor.media.substrato
ss.intra=sum((tamboril$altura-tamboril$media.grupo)^2)

#somatrio do desvio quadrtico entregrupo
media.substrato
media.altura
ss.entre=6*sum((media.substrato-media.altura)^2)

#Calculo do desvio mdio intra grupos dividindo o ss.intra 
#pelo gl intra, e o desio mdio entre grupos dividindo o ss.entre 
#pelo gl entre.
gl.total=60-1
gl.entre=10-1
gl.intra=gl.total-gl.entre
dm.entre=ss.entre/gl.entre
dm.intra=ss.intra/gl.intra
#Calculo da razo de variancia F
razao=dm.entre/dm.intra
#Probabilidade da diferena
prob=pf(razao,9,50, lower.tail=FALSE)
prob
porc=100*(ss.entre/ss.total)

####Exerccio 6.3#############

mudas<-read.table("altura-mudas.csv", head=T, sep=",")
head(mudas)
tamboril<-mudas[mudas$especie=="tamboril",]
tamboril$substrato<-as.factor(tamboril$substrato)
str(tamboril)
aov.tamboril<-aov(tamboril$altura~tamboril$substrato)
aov.tamboril
summary(aov.tamboril)

